Magnetically-induced reconstructions of the ground state in a few-electron Si 

quantum dot 
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We report unexpected fluctuations in the positions of Coulomb blockade peaks at high magnetic 
fields in a small Si quantum dot. The fluctuations have a distinctive saw-tooth pattern: as a function 
of magnetic field, linear shifts of peak positions are compensated by abrupt jumps in the opposite 
direction. The linear shifts have large slopes, suggesting formation of the ground state with a non- 
zero angular momentum. The value of the momentum is found to be well defined, despite the 
absence of the rotational symmetry in the dot. 

PACS numbers; 73.23.Hk, 85.30.Wx, 85.30.Vw, 85.30.Tv, 71.70.Ej 



The basics of Coulomb blockade (CB) phenomena can 
be understood within the so-called "orthodox theory" 
iQl . In this theory, electron-electron interactions are hid- 
den in the charging energy and the electrostatic cou- 
pling is assumed to be independent ofO the nature of 
the ground state, namely on the particular distribution 
of electrons inside the dot. Early experiments on large 
vertical quantum dots have already revealed significant 
deviations from the " orthodox theory" , noting that sev- 
eral electrons can enter the dot at almost the same energy 
1^. The deviations were later attributed to localization 
of electrons in local minima of the confining potential . 
In a smooth confining potential, magnetic field forces a 
redistribution of charges within the dot to form quantum 
Hall edge states [Q. Charge redistribution is a focus of 
much theoretical work, especially in the regime of high 
magnetic field Q . A competition between the attractive 
confining potential and the repulsive electron-electron in- 
teractions is expected to produce a rich variety of exotic 
patterns of charge distribution In a recent experi- 
ment, high field (filling factor v < I) redistribution of 
charges in a large vertical dot > 0.5 /Ltm has been re- 
ported 0. 

So far, most of the experiments were performed on two- 
dimensional dots with weak and smooth confining poten- 
tial electrostatically created by gating. In such dots, the 
magnetic field B dependence of energy levels is domi- 
nated by orbital effects even at low B. Recently, it has 
become possible to investigate quantum dots in a differ- 
ent regime of very small size, strong confinement, and 
strong local disorder . This regime is realized in three- 
dimensional Si dots with confining potential provided by 
a sharp Si/Si02 interface. In these dots states with dif- 
ferent angular momenta are mixed due to the absence 
of rotational symmetry, eliminating the linear in B term 
of the orbital energy, and the parabolic B^ term is sup- 
pressed by strong confinement. Thus, the B dependence 
of the energy levels is expected to be simple and to consist 
of linear Zeeman shifts. Indeed, earlier we demonstrated 
that the shifts of CB peak positions are dominated by 
the Zeeman effect [01. In this work we focus on small 



fluctuations in the peak position, which appears at high 
B and low T in a Si quantum dot with a few electrons, 
N K A — 6. The fluctuations have a distinct saw-tooth 
pattern: as a function of B, fast linear shifts in the peak 
position are interrupted by abrupt jumps in the opposite 
direction. We show that the linear segments of the fluctu- 
ations have orbital origin, indicating formation of states 
with a relatively well defined non-zero angular momen- 
tum close to 1 or 2. This is an unexpected result, taking 
into account the geometry of the sample. Also, we show 
that the jumps arc intrinsic to the dot, presumably due 
to some magnetically induced rearrangement of charges 
inside the dot. We discuss the data within a simplified 
single-particle picture. However, the observed phenom- 
ena is clearly a many-body effect, and a proper treat- 
ment of interactions, strong confinement and disorder is 
needed. 

The small Si quantum dot is lithographically defined 
from a silicon-on-insulator wafer, and a sample layout is 
shown schematically in the inset in Fig. |l| A detailed de- 
scription of sample preparation can be found in Ref . . 
The dot is three-dimensional, asymmetric, and is elon- 
gated in the current flow direction. A poly-Si gate is 
wrapped around the dot and is used to control the num- 
ber of electrons in the dot starting from iV = 0; the gate is 
separated from the dot by 500A of Si02. The dot is con- 
nected to two two-dimensional source and drain contacts 
via tunneling barriers; the coupling is weak and even at 
the lowest temperature, T — 60 mK, the CB peaks re- 
main thermally broadened. From excitation spectra we 
deduce single-particle energies to be 1 — 5 meV, compa- 
rable to the charging energy 5 — 7 meV. 

The uniqueness of our dot is its small size 100 — 200 
A long and < 100 A in cross section) and the extremely 
strong confining potential provided by the Si/Si02 inter- 
face (~ 3 eV). The strength of the confinement is clearly 
seen from the analysis of the first peak, see Fig. 2 of 
Ref. ||l^. The first energy level has a weak parabolic 
B dependence due to magnetic confinement {hujc)'^ / Eq, 
where u)c = eB/m* is the cyclotron frequency. Char- 
acteristic energy Eq depends on the direction of B: 
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FIG. 1. (a) position and (b) amplitude of peak 5 plotted as 
a function of magnetic field for two field directions: in-plane 
field _B|| along the current direction (solid line) and B± per- 
pendicular to the Si wafer (dots). Inset shows a schematic of 
the sample. The data was measured at T = 60 mK with 10 
fiY ac bias. Dashed lines have slopes ^g*fiB, g* = 2. 

Eq k, 100 mcV for B applied perpendicular to the sam- 
ple, B±, and there is no measurable shift of the first level 
for in- plane B aligned with the current direction, i?| | . For 
both B directions the magnetically induced confinement 
is smaller than the Zeeman energy in the experimental 
range of B. 

The number of electrons in the dot can be tuned be- 
tween and 30 with the gate voltage Vg. The overall 
B dependence of the energy levels in this sample has 
been analyzed previously in Ref. Evolution of the 
fifth CB peak with B is plotted in Fig. |l|a. According 
to the previous analysis, the V^{B) curve consists of 3 
segments, each having k, ^g*iiB/a slope, where g* = 2 
is the ^-factor in Si, = eh/2mQ is the Bohr magneton 
and a ^ 14 mV/meV. The two kinks at B « 2.5 T and 
B Ri AT have been identified as crossings of Zeeman split 
spin-up and spin-down levels, originating from different 
single-particle energy states. Indeed, average peak posi- 
tion does not depend on the direction of B. However, 
there is a noticeable difference between the V^{B^) and 
VP{B±) curves at high B: changes linearly with B^ 
but fluctuates as a function of B±. The appearance of 
these fluctuations is accompanied by a dramatic increase 
of the peak amplitude by more than two orders of 
magnitude, as shown in Fig. |l|b. The fluctuations are 
most pronounced at our base temperature of 60 mK and 
are still observable at 200 mK, while at T « 1 K they are 
washed out completely. No fluctuations were observed 
for > 20 or < 4. The fluctuations are highly repro- 
ducible: the measured VJ and from different scans 
are identical if the scans are performed within the same 
cooldown; patterns from different cooldowns are similar, 
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FIG. 2. Evolution of peaks 4, 5 and 6 as a function of 
perpendicular magnetic field On the top panels, a series 
of G vs. Vg curves are plotted, each linearly offset with B^. 
The bars are 4 scales. In the bottom panels, peak positions 
are extracted and plotted as a function of B\_ for the same 
range of Bx as in the corresponding top panels. Solid and 
dashed lines have the slopes ^Uloc/B and htOc/B. 

although positions of the jumps may vary. 

Fluctuations of VJ for peaks 4, 5 and 6 are magnified 
in Fig. ||. The fluctuations have a distinctive saw-tooth 
pattern: increases linearly with B^, then decreases 
abruptly. Most of the linear segments have slopes close 
to 4 mV/T (0.3 meV/T) or 8 mV/T (0.6 meV/T). These 
slopes are much larger than ^g*fiB — 0.06 meV/T and 
are grouped within 15% around ^huJc/B and Tiujc/ B. The 
linear shifts are interrupted by abrupt jumps of the peak 
position, which occur within < 20 mT. These jumps com- 
pensate the linear shifts and, on average, the peaks follow 
the weak vs B\\ field dependence due to the Zeeman 
shift of energy levels. 

Orbital effects can be distinguished from spin effects by 
studying the evolution of peak positions in a tilted mag- 
netic field. In Fig. |^, VJ of peak 5 is plotted as a function 
of i3sin(8) for different angles 8 between B and /. Lin- 
ear peak shifts depend on the perpendicular component 
B±^ « i3sin(0), rather than on the total filed B (solid 
lines provide a guide for the eye). Curves at different 
angles do not scale with sin(8) exactly, as is expected 
for a three-dimensional structure. Taking into account 
that the energy shift due to the magnetic confinement 
is negligible, we attribute the slopes « ^Tiloc/B to the 
formation of states with a non-zero angular momentum 
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FIG. 3. Position of peak 5 as a function of B± — B ■ sinO 
plotted for different angles Q. The angle is defined in the 
inset. The lines have a slope ^TiuJc/B and are guide to the 
eye. 

m 1. Segments with larger slope « huJc/B have also 
been observed and we attribute them to the formation of 
states with m « 2. 

Fast linear shifts of as a function of B± are inter- 
rupted by abrupt jumps in the opposite direction. These 
jumps occur within AB < 20 mT and, if attributed to 
a state with large angular momentum, would correspond 
to an unreahstic m > 10. Ground state (GS) energy of 
the dot should be a continuous function of all variables, 
including B, unless the GS is bi-stable. Bi-stability of the 
GS should reveal itself through a hysteresis of the conduc- 
tance in Vg or B scans. Experimentally, the traces were 
found to be identical, independent of the scan direction. 
Thus, jumps in the peak position should reflect either 
abrupt changes in the environment or abrupt changes in 
the electrostatic coupling between the dot and the envi- 
ronment. 

In the CB regime, conductance is non-zero only when 
the electrochemical potential of the dot equals the elec- 
trochemical potential of the leads. Peak position is 
determined from the condition Ifllll 



e[(jV - 1/2) - CgVP] 

Cy^ 



+ [E{N)^E{N-l)] + e^, 



where e^/2Cs is the Coulomb energy for one additional 
electron, Cg and Cs denote respectively the gate and the 
total capacitances, E{N) is the kinetic energy of the state 
with N electrons, Ep is the Fermi energy in the leads, and 
(j)ex is a potential induced by external charges. According 
to the above equation, there are three possibilities for a 
jump of V^: i) an abrupt change in Ep, ii) a discrete 
change of the background charge distribution and, thus, 
4>ex 1 or iii) an abrupt change in the electrostatic coupling 
Cg. It is easy to rule out Ep as a source of the jumps, 
since depopulation of Landau levels should lead to the 
jumps of in the opposite direction. 
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FIG. 4. Gray-scale plots show conductance as a function of 
Vg and _Bj_ near peak 5. The data was taken by scanning Vg 
at different fixed values of Bx. In (a) the amplitude of each 
Vg jump is extracted. In (b) conductance as a function of Bx 
was measured at a fixed Vg — 0.7005 V, marked by the dashed 
horizontal line in the gray-scale plot. Note the logarithmic 
scale of G. The vertical dotted lines separate ground states 
with different number of electrons along the curve. The data 
in (a) and (b) were collected during different cooldowns. 



It is appealing to attribute the jumps to a magnetically 
induced depopulation of a charge trap. However, such 
an explanation is inconsistent with the data. One has to 
assume the existence of a large two-dimensional trap ca- 
pacitively coupled to the dot, which does not participate 
in the charge transport. Qualitatively, there are rather 
general arguments against such a scenario. If jumps are 
related to the magnetic confinement of electrons in the 
trap, we expect the frequency of the jumps to increase 
with B, and there should be no low-field cut-off for their 
appearance. Experimentally, the jumps appear suddenly 
at B 5 T and the frequency of the jumps slightly de- 
creases with B; there are no jumps at lower fields. In Si 
nanostructures, traps are formed inside the oxide layer 
|]l2| , and it is hard to explain their existence only at 
high B± > 5 T and in a limited range of gate voltages 
0.5 Y < Vg < 1 V. Quantitatively, in order to account 
for the data, one has to assume the existence of a) sev- 
eral traps, b) all within the tunneling range from the 
dot, inconsistent with the geometry of the sample. Elec- 
trostatically, discharging of a trap leads to a fixed shift 
AVP = —Cc/C^ ■ e/Cg, where Cc is the dot-trap cross 
capacitance, and Cs and are the corresponding to- 
tal capacitances. In Fig. ^a, values of the shifts are 
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extracted for each jump in the range 5 T < < 8.2 T. 
First, there is a positive shift at B± = 7.4 T, inconsistent 
with the trap discharging model. Second, amphtude of 
the negative shifts ranges from 0.2 to 1.0 mV with no sys- 
tematic pattern. Thus, a single trap cannot account for 
the observed oscillations. Subsequent discharging of sev- 
eral traps can be ruled out using the following arguments. 
In the Vg — B plane the number of electrons in the dot 
differs by one across the jump, as shown schematically in 
the inset in Fig. Following the analysis of a two-dot 
system Q, there should be a peak in the conductance 
along the jump, unless the trap and the dot are connected 
by tunneling (in this case an electron tunnels between the 
dot and the trap and the total number of the electrons in 
the dot-trap system remains constant). Experimentally, 
there is no peak in conductance if B is swept across the 
jump at a fixed Vg, as shown in Fig. Qd. An assumption 
that several large traps are in a close proximity of the 
dot is inconsistent with the small size of the device. 

It is very suggestive that all the observed effects, 
namely jumps in peak position, enhanced conductivity, 
appearance of large slopes and an apparent mutual com- 
pensation of shifts and jumps, are related and originate 
from the dot. Presumably, the jumps are related to the 
magnetically induced rearrangement of electron density 
inside the dot. In our dots, we expect strong fluctua- 
tions of the confining potential, because Si/Si02 interface 
roughness directly translates into large fluctuations of the 
local potential. At low B fields the electron wavefunction 
is spread over the entire dot. At high B fields the extent 
of the wavefunction is determined by the magnetic length 
Im = {h/eBY^^. When Im becomes smaller than the 
average distance between electrons. Coulomb repulsion 
favors localization of electrons. We expect fluctuations 
of the local potential to facilitate the redistribution of 
charge. Strong increase in the peak height hints that elec- 
trons are rearranged closer to the dot boundary, where 
their wavefunctions have larger overlap with electrons in 
the leads. It has been pointed out in Ref. jlj] that elec- 
trostatic coupling to the gate Cg is not necessarily a con- 
stant, but depends on the particular distribution of elec- 
tron density inside the dot. In order for a peak to shift by 
1 mV, only a small change ACg/Cg « AI///V7 < 0.2% 
is required. 

The most unexpected and surprising result is the ap- 
pearance of states with non-zero angular momentum. In 
the absence of rotational symmetry, states with differ- 
ent angular momenta are mixed, and angular momen- 
tum is not a good quantum number. However, the dot 
asymmetry becomes less important if an electron can 
complete a classical cyclotron orbit faster than the time 
r = L/ y^2E/m* needed to traverse the dot. For an 
electron with kinetic energy E — A meV in a dot of size 
L = 200 A, the cross over ^ < r should occur at i? « 5 
T. Note that E k, N/S. increases with the number of elec- 
trons (here A is the level spacing) , leading to the increase 



of the cross over field with N. We speculate that, for a 
few electrons in the dot, high B provides a mechanism 
to suppresses the mixing of different angular momentum 
states. 

Electron-electron interactions in our sample are rather 
strong {e^ /A-KetQT = 12 meV for r = 100 A). They reveal 
themselves in the spontaneous polarization of the N = 6 
ground state at low fields and in the suppression of the 
corresponding CB peak due to spin blockade ^. Thus, 
a many-body description of the ground state is required. 

To summarize, we explored a quantum dot in a new 
regime of small size, strong confinement and strong 
electron-electron interactions. We observed fluctuations 
in the CB peak positions at high B for a few-electron 
states. The fine structure consists of fast linear shifts, 
followed by abrupt jumps in peak positions. We argue 
that the fluctuations are intrinsic to the dot and reflect 
changes in the many-body wavefunction. We attribute 
jumps to the magnetically induced spatial rearrangement 
of charges inside the dot. The linear segments have large 
slopes, reflecting unexpected formation of states with rel- 
atively well defined non-zero angular momentum. The 
observed phenomena is a many-body effect and is clearly 
beyond the description of a single-particle picture. 
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reading of the manuscript and valuable comments. The 
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